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Abstract 

A 2D model describing depinning of an interface from a rough, self-affine 
substrate, is studied by transfer matrix methods. The phase diagram is deter- 
mined for several values of the roughness exponent, £5, of the attractive wall. 
For all Cs > the following scenario is observed. In first place, in contrast to 
the case of a flat wall (£s = 0), for wall attraction energies between zero and 
a Cg-dependent positive value, the substrate is always wet. Furthermore, in 
a small range of attraction energies, a dewetting transition first occurs as T 
increases, followed by a wetting one. This unusual reentrance phenomenon 
seems to be a peculiar feature of self-affine roughness, and does not occur, e. 

g., for periodically corrugated substrates. 
PACS number (s): 68.45.Gd 
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I. INTRODUCTION 



Wetting phenomena occur when, e. g., a layer of liquid phase coexisting with its vapor 
grows macroscopically over an attractive solid substrate Q. While wetting in pure systems 
is relatively well understood by now, the effects of disorder in interfacial phenomena @-[5| 
pose many challenging issues and are the object of active research. Of particular interest 
is the effect of geometric surface disorder (roughness) on the location and the nature of the 
wetting transition. An extensively studied type of roughness, also in view of its experimental 
realizability B , is that occurring when the average height fluctuation in a sample of longi- 
tudinal linear size L, wl, scales like wl ~ L^ s . This self-affine scaling disorder is globally 
characterized by the roughness exponent Cs> which can be expected as the only relevant 
substrate parameter possibly affecting universal features of the wetting transition. 

The relevance of Cs f° r both complete and critical wetting transitions has been extensively 
studied in the last years As far as critical wetting is concerned, a recent study by the 
present authors || has put in evidence the fact that self-affme roughness can produce a 
change from continuous to first-order transition in systems with short-range forces, when 
the roughness of the substrate, Css exceeds the intrinsic roughness of the interface in the 
bulk, Co- On the other hand, for Cs < Co 5 no change is expected in the nature and universality 
class of the continuous wetting transition, with respect to that on flat substrate. 

The possibility that substrate roughness drives a continuous wetting transition first- 
order has been subsequently discussed in the context of Landau-type mean field approaches 

A possible modification of the nature of the wetting transition is not the only effect of 
self-affine substrate roughness on wetting. As we show here, roughness produces modifica- 
tions of the phase diagram of the interface, which can be quite dramatic and important in 
experiments and applications. 

In the present paper we study systematically the wetting phase diagram of a gener- 
alization [|9] of the Chui- Weeks model [DJ with a rough attractive boundary in 2D. By 
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this study we produce evidence of some remarkable and definitely unusual features of the 
phase diagram, which, at a qualitative level, should be considered as generic for wetting on 
self-affme rough substrates. The most notable feature of the phase diagram is a reentrance 
phenomenon for wetting in the whole range of roughnesses (0 < (s < !)■ This reentrance, 
which implies a dewetting followed by a wetting transition as the temperature is raised, 
occurs both in regimes when the transitions (both dewetting and wetting) are critical, and 
when they are first-order. 
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FIG. 1. Example of rough substrate wall (continuous path) and interface configuration 
(dotted path). 

This paper is organized as follows. In the next section we introduce the model and 
describe the transfer matrix method we apply. In section 3 we dicuss how the phase diagram 
is calculated and illustrate in detail the main results. The last section is devoted to further 
general considerations and to conclusions. 

II. MODEL AND TRANSFER MATRIX 

We consider here the same generalization of the Chui-Weeks interfacial model discussed 
in ref. ||. Let us denote by x and y the integer coordinates of points on a square lattice. 
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The interface configurations are self-avoiding paths (partially) directed in the x direction 
(see Fig. 1). Each interface step parallel to the x axis is located by giving the ordinate 
y = h x of its left-hand extremity. In this way the configuration is determined by giving h x , 
Vzc. We impose the following extra restrictions on the set {h x }. First of all h x > S x where 
{S x } represents the directed profile of the substrate wall. {S x } are sets of integer ordinates 
determining the wall configuration with the same conventions applied to {h x }. An extra 
constraint on possible interface configurations is given by h x+ \ — h x = 0,±1, while for S x 
we choose S x+ \ — S x = ±1. We impose such constraints for computational convenience: 
their removal or modification would not change the main qualitative features of the phase 
diagram. 

The sets {-S^} are generated by a random sampling procedure inspired by Mandelbrot 



12|| and described in Appendix. This procedure produces directed paths in 2D obeying the 
restrictions described above and satisfying the scaling relation 

\S X+ Ax - S x \ ~ |Ax| Cs (1) 

In the last equation the bar indicates statistical average with respect to the above men- 
tioned sample of {S x }. Averages over wall configurations, i. e. over different {S x }, will be 
considered quenched and denoted by overbars below. 

The interface Hamiltonian is of the form 

x 

H = X^K 1 + \ z x ~ Zx-i + s x \) - ud Zxfi ] (2) 

with z x = h x — S x and s x = S x — S x _i. In eq.(||) e (e > 0) is the energy cost of any interface 
step and — u (u > 0) is the energy gain of an interface contact with the attracting wall. The 
sum in (0) is performed up to X, which represents the length of the interface projection on 
the x axis. It should be noted that in our model only the horizontal steps of the interface 
paths in contact with the wall are prized by energies —u. This choice is again not mandatory. 
Different conventions would not change the basic qualitative results. 

At a finite temperature T the fugacities to = e _£//T and k = e u ^ T are associated with 
each (horizontal or vertical) step of the path, and to each horizontal step on the wall, 



respectively. Consequently, given a wall profile, the partition function associated to all the 
interface configurations covering a distance X can be written in the form 

Zx _ e ~T[ e E» ( 1 +l*-««-i+*-l)- u E**-.o] = u x to n± k nc (3) 

{21} {z x } 

where the sum is done over the ensemble of all the directed paths (determined here by 
{z x }) compatible with the choosen wall. n± and n c are the number of vertical steps of the 
interfacial path and the number of horizontal steps on the wall, respectively. Note that the 
total length of a path is given by L = X + n±. 

The partition function ([3D can be more conveniently expressed in the form 

2*=^£(n T ..) (4) 

l,z \x=l ) 1 >Z 

where the T Sx are transfer matrices defined as 

(T s ) m>n [<^m,n— s ~t~ ^i^m.n— s— 1 $m,n— s+l)]^ (5) 

with m,n > 0. The function O m ec i- ® can be used to enforce particular initial conditions 
for the interface paths; if we choose paths with an extremity on the wall, we put 4>o(z) = 5 Z) o- 
With the above definitions, a wall profile corresponds to a particular sequence of factors 
T Sx in the product of eq.(|j). Correspondingly, for asymptotically large systems, the partition 
function Z x can be expressed in terms of the largest Lyapunov eigenvalue X max as Z x ~ 
{w\ ma x) X , for X — > oo, where Jl3| 



^max — lim 

A^OO 



X 



;n T ^)^n/ii0oi 



(6) 



x=l 

with||0||=E 2 0W-[] 

For walls with some periodic geometry with period X p the matrix product in (0) can be 
written as [T Xp } x ^ Xp (X assumed an integer multiple of X p ) with T Xp = n^=iT S:E . In this 



1 Note that our definition of norm is legitimate because 4>o(z) and the elements of the T Sx 's matrices 
are > 0. 
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case the calculation of \ ma x is equivalent to the calculation of the largest eigenvalue, Ax p , 
of the matrix T Xp . Consequently, X m ax = [^■x p ) 1 ^ Xp - 

In numerical calculations it is useful to introduce the normalized vectors ip x defined by 
the recursion relation 

f x = —T S J X ^ (7) 

n x 

with n x = \\T 8x iJ) x -i\\ and ipo = <fio- It is straightforward to see that the z-th component 
of the vector ip x corresponds to the probability that the path at a; is at a distance z from 
the wall @||. This suggests our definition of H^ll- F° r a given wall profile, the above 
definitions allow to express the Lyapunov eigenvalue (0) as 



^max — lim 

X — >oo 



X 
Lx=l 



1 X 



x=l 



The quenched, dimensionless free energy density given by limx-»oo — hi Zx/X can be written 
in the following form 



1 x 

f = -lnuj\ max = -hiu - lim — V lnn x . (9) 

X^ooX x=1 

If X is chosen large enough, the average over quenched wall disorder for the second term on 
the right hand side requires only a rather limited sample of wall configurations. This is due 
to a self-averaging property of /, which clearly manifests itself in the numerical results. 



III. DEPINNING TRANSITION AND PHASE DIAGRAM 

In the discussion below we implicitly make u and T dimensionless by dividing them by 
e. The wetting transition occurs because, e. g., at a given T, the interface can be bound 
to the wall of the substrate only for sufficiently high values of u. In the case of ordered flat 
walls, i. e. {S x = constant Va;}, the value of u above which the interface is pinned can be 



easily calculated [14,11,21 to be 



u c {T) = Tin [(1 + 2 exp(-l/T))/(l + exp(-l/T))] (10) 



with u c (0) = 0. We denote by Pq the average fraction of interface horizontal steps on the 
wall, 

P = Km (n c )/X, (11) 

X — >oo 

with brackets indicating canonical thermal average. One can check that Pq vanishes contin- 
uously and linearly when the line u = u c (T) is approached from above. 

When dealing with our random substrates the calculation of Pq or / for each {S x } can not 
be done exactly in a semi-infinite geometry, and truncations must be made. To minimize 
effects of the finite size of the transfer matrices, in our calculations we only considered 
matrix sizes much larger than the mean square perpendicular width of the self-affme walls. 
This width can be defined as AS ± = (S* - S 2 ) 1 ' 2 ~ with S 2 = (l/X)Y JX S 2 x and 

S = (1/X) J2 X &x- m practice we used transfer matrices as large as 10 4 x 10 4 in the roughest 
case, corresponding to (s — l n 12/ In 32 ~ 0.717. With this roughness X = 10 5 was reached. 

In addition one has to average over different {S x } in order to get Pq and /. In practice 
we could sample at most 10 or 15 different {S x } in the most favorable cases, due to the large 
X's needed to extract precisely P . Fortunately, as mentioned above, our large X values, of 
the order of 10 5 , lead to a very high degree of self-averaging in quantities like / and Pq. 

At fixed T, as the transition is approached from above (viz. u > u c (T)) the interface free 
energy density @ is negative and increasing, with decreasing u; at u = u c (T) it matches 
the bulk interface free energy density fb u ik = — lnu;(l + 2u). Thus, it is useful to consider 
the interface excess free energy density, Af = f — fbuik, which, by means of eq. @, can be 
expressed as 

Af = -hiX^aZ + ln(l + 2a>) = -InA^ (12) 

This last equation defines X' max as the Lyapunov eigenvalue of transfer matrices defined as 
T' s = [1/(1 + 2c<j)]T s , which were directly used in our calculations. 

Thus the position of the depinning transition corresponds to the vanishing of Af. Also 
the interface contact probability can be evaluated using X' max , taking into account that 
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FIG. 2. Interface phase diagram for rough self-afhne walls in the T-u plane (e = 1 
is assumed). The curves u = u c (T) are shown for five different values of the roughness 
exponent (s- The light continuous line corresponds to u = u c (T) for a flat wall as given by 
eq. (0). 



The calculation of Pq offers an alternative way of locating dewetting or wetting transi- 
tions, by identifying the conditions under which the quantity first becomes zero. A numerical 
study of how P$ approaches zero can also provide informations on whether these transitions 
are continuous or discontinuous ||. We used this criterion, together with an analysis of the 
way in which Af vanishes to clarify the first or second order nature of the transitions. 
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Fig. 2 summarizes the results we have obtained by a systematic calculation of / as 
function of T and u, for 5 different values of the (s- The curves in the figure represent 
the behavior of u c versus T. u c was determined numerically as the value of u below which 
< 0.0001 as u was changed in steps of 0.001. By using the definition of uo and k the 
u = u c (T) curves can be converted into the curves k = k c {uj). Note that u c (0) ^ implies a 
power-law divergence of k c (u>) for uj — > 0: k c ~ o; _Uc ( ). 

Looking at the curves for rough walls in Fig. 2, we note two main differences from the 
flat case: 

a) It is not possible to pin an interface to a rough wall with a vanishing contact energy. 
As T — > the minimal contact energy to pin an interface, u c (0), is finite and increases 
as the wall roughness, i. e. £g, increases. 

b) For each ( s , there is a temperature, T R (( S ), below which u c is a decreasing function 
of T. Note also that for all ( s < 1 u c (T) approaches u c (0) with zero slope. 

A surprising consequence of b) is that if, as in an experiment, we monitor interface 
behavior at fixed u, by varying T, we have to distinguish among three kinds of regimes: 

I) For u < u c (T R ) interface pinning is impossible, no matter how low T is. The substrate 
is wet at all T, and no transitions take place. 

II) for u > u c (0) as the temperature is increased the interface passes from a pinned to a 
depinned state at some Tw- Thus, the substrate is partially wet for T < Tw and wet 
for T > T w , as for smooth walls. 

Ill) for u c {T R ) < u < u c (0), as the temperature is increased, the interface undergoes 
two transitions: 1) at some temperature T D < T R we find an unexpected dewetting: 
the interface which is depinned for T < Tp becomes pinned at T = T^; 2) at some 
Tw > T R a more usual depinning transition follows. Thus, the substrate is wet at very 
low T, then it dewets, and eventually it wets again at high T. 
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Unfortunately, a detailed study of the behavior of Tr for (s approaching zero is not 
feasible due to the necessity of generating extremely long walls in order to distinguish, e. g., 
(s = 0.1 from (s = 0. However, our results suggest rather clearly that Tr approaches zero 
for both (s — > and — > 1. Thus, in these two limits the reentrance desappears. 

Another interesting aspect of the phase diagram is that connected to the nature of the 
transitions involved. The continuous or discontinuous character of the wetting transitions 
upon varying was already discussed in ref. |§ by analyzing the way in which Pq approaches 
zero for u — > u c (T). There we found that when (s ~ 1/2 the transition remains second- 
order and most likely belongs to the same universality class as the flat case, i. e. Pq ~ 
(u — u c )^, with ip = 1. When (s exceeds 1/2 there is clear evidence that depinning occurs 
disco ntinuously. Co — 1/2 is the roughness exponent of the interface in the bulk and it makes 
sense that this is the precise border value of separating the two regimes. 

A natural question is whether, in the range u c {Tr) < u < m c (0), the dewetting transition 
occurring at low temperature is of the same kind as its wetting counterpart at higher T. 
Following the lines of ref. || we made a systematic study of the dewetting transition for 
two (s values, respectively below and above ^ = 1/2. In the first case ((s = 0.4) we found 
evidence of a continuous dewetting, while in the latter (( s = 0.6) it appeared discontinuous. 
Thus, in spite of the fact that the dewetting transition occurs only for some u ranges, and 
at much lower temperatures, it seems that its character could be the same as that of the 
corresponding wetting transition. Of course, a more precise definition of the limits within 
which the above identity of transition orders applies would require an extensive systematic 
exploration, which is beyond the scope of the present work. 

IV. CONCLUSIONS 

The results presented in the previous section are somewhat unexpected and are worth 
discussing further. 



10 



< > 

Xp 



FIG. 3. Example of a periodically corrugated wall on a square lattice (heavy line). The 
wall period is X p = 10. The light straight horizontal line corresponds to the bound interface 
ground state at T = and u > 0. 

The first important fact is that u c (0) > for ( s > 0. At T = 0, in order to decide whether 
the interface is pinned or not, we need only to compare the ground state energy in the bulk 
with the lowest energy of a state in which the interface is bound to the substrate. In the 
bulk the state of lowest possible energy is clearly given by a straight configuration (n_i_ = 0). 
A bound state will have an energy relative to this unbound ground state equal to n_i_ — un c 
(assuming e = 1). n_i_ and n c of course depend on the wall configuration to which this bound 
state refers. Clearly u c (0) is determined by the condition under which this energy difference 
between the two states vanishes: u c (0) = limx^oo n±/n c . The fact that bound ground state 
configurations satisfy this limit condition with u c (0) > is a nontrivial property of self-affine 
substrates. On a periodically corrugated substrate with average horizontal orientation (like 
that sketched in Fig. 3) this limit property would not be satisfied. In that case, for u very 
close to zero, the bound ground state configuration is one in which n± = 0, corresponding 
to a straight interface touching the attractive tips of the periodically corrugated wall (see 
Fig. 3). Thus, we would have n± = and n c ^ 0, and, consequently, u c (0) = 0, like in a 
flat case. We conclude that a remarkable property of self-affine substrates is that they can 
support ground state interface configurations with limx->oo n±/n c > 0. 

Notice further that for a periodically corrugated substrate the reentrance phenomenon 
with dewetting preceeding wetting is not possible, since in that case the curve u = u c (T) 
can not be decreasing in the neighborhood of T = 0. We verified by explicit calculations 
that for a wall as in Fig. 3, u c (T) is in fact never decreasing on the whole T axis. 

An even more remarkable property of the self-affine substrate, for which u c (0) > is 
clearly a necessary but not sufficient condition, is the monotonically decreasing character of 
the curve u = u c (T) in the interval (0, T R ). This feature implies that, as soon as T > 0, an 
interface can be more easily bound to the rough substrate. This clearly shows that there 
is a very nontrivial energy-entropy interplay in the pinning mechanism when self-affine 
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roughness is involved. 

We also compared the nature of dewetting transitions to the corresponding high-T wet- 
ting ones. We got only preliminary results on this issue. These results suggest the possibility 
that, once a first-order character prevails for wetting, the same applies to dewetting as well. 
This would mean that geometry alone is the crucial factor in determining the nature of 
transitions on rough substrates. 

The model calculations we presented here are of course limited to 2D and to strictly 
short-range forces. An extension in 3D is computationally unfeasible, and also the inclusion 
of long-range potentials would pose serious additional difficulties in our calculations. In 3D 
we expect that the main features of the phase diagram would persist. Of course a major 
difference in 3D would be the character of the transitions. In 3D the interface has roughness 
Co = in the bulk [0 . Thus, following the conclusions of ref . || , we should probably expect 
first-order dewetting and wetting transitions for all (s > 0. 

Concerning the effect of long-range forces, which should certainly be included in more 
realistic calculations to compare with experiments, we can only conjecture that they would 
not modify the main result obtained here, i. e. the reentrance. However, there is at least 
one instance, that of interfaces in superconductors |TJ|, in which a short range description 
is fully adequate. 
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APPENDIX: GENERATING RANDOM SELF AFFINE PATHS IN 2D 

The procedure we describe here is a random version of a deterministic algorithm by 
Mandelbrot ||12||. Given two even integers, p and q, with p < q, the procedure allows us to 
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construct iteratively a partially directed path with a roughness exponent (s = hip/ lng and 
with X = q n after n iterations. 

First we consider a set {a} of vectors with q components. For each vector a in the set 
(q + p)/2 components, chosen at random, are set equal to 1, while the remaining (q — p)/2 
are put equal to —1. Once the set {a} contains a sufficiently large number of such elements, 
the construction of a wall profile proceeds iteratively. In the first iteration we consider a 

— * 

vector Pi with q components and, after choosing an element a at random in the set {a}, we 

— * 

set Pi = a. 

In the second iteration we construct a vector p 2 with q 2 components. Once chosen q 
vectors dj (j — 1, 2, . . . , q) at random in {a}, the components of p 2 are determined according 
to the rule: 

/M0'-i)g + i) = A0>;W ( A1 ) 

with i — 1, 2, . . . , q and atj(i) indicating the i-th component of the vector dj. 

The last equation can obviously be iterated for the construction of vectors p n with q n 
components: 

Pn((j ~ I)?"" 1 + l) = Pn-lUhjii) (A2) 

with j — 1,2, ... , q n ~ l and % = 1, 2, . . . , q. 

At any iteration n a directed path defined by the sequence of integers S x , as explained 
in section 2, can be obtained from p n by 

S x+1 = S x + p n {x) (A3) 

where we usually set Si = 0. Larger n implies a large horizontal size X = q n of the path. The 
accuracy of the self-affine average scaling of the paths, as detectable on the basis of eq.(l), 
increases with increasing n and the number of walls considered. For q ranging between 8 
and 32 we have seen that when the longitudinal length of the paths X is of the order of 10 5 
the computed roughness exponents of a single generated path coincides with the theoretical 
(s = lnp/ lng within 1%. 
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